Decay laws for three-dimensional magnetohydrodynamic turbulence 
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Decay laws for three-dimensional incompressible magnetohydrodynamic turbulence are obtained 
from high-resolution numerical simulations using up to 512 3 modes. For the typical case of finite 
magnetic helicity H the energy decay is found to be governed by the conservation of H and the 
decay of the energy ratio V = E V /E M . One finds the relation (E 5/2 /eH)Y 1/2 /(l + T) 3/2 = const, 
e = —dE/dt. Use of the observation that F(t) oc E(t) results in the asymptotic law E ~ t~ ' 5 in 
good agreement with the numerical behavior. For the special case H — the energy decreases more 
rapidly E ~ t~ , where the transition to the finite- H behavior occurs at relatively small values. 

PACS: 47.65+a; 47.27Gs; 47.27Eq 



Many plasmas, especially in astrophysics, are charac- 
terized by turbulent magnetic fields, the best-known and 
most readily observable example being the solar wind. 
The convenient framework to describe such turbulence 
is magnetohydrodynamics (MHD). Here one ignores the 
actual complicated dissipation processes, which occur on 
the smallest scales and would usually require a kinetic 
treatment, assuming that the main turbulent scales are 
essentially independent thereof. Instead dissipation is 
modeled by simple diffusion terms. If, moreover, interest 
is focussed on the intrinsic turbulence dynamics, one can 
also ignore the largest scales in the system, which depend 
on the specific way of turbulence generation, restricting 
consideration to a small open homogeneous domain of 
the globally inhomogeneous turbulence. 

Homogeneous MHD turbulence has become a 
paradigm in fundamental turbulence research, which has 
been receiving considerable attention. It is well known 
that 2D and 3D MHD turbulence have many features 
in common concerning, in particular, the cascade prop- 
erties. In both cases there are three quadratic ideal 
invariants: the energy E = | J (v 2 + B 2 )dV, the cross 
helicity K = J v • Bcfl/', and a purely magnetic quantity, 
the magnetic helicity H = J A • TSdV in 3D and the 
mean-square mag netic potential = J ^ 2 dV in 2D, 
which both exhibit an inverse cascade. Many theoretical 
predictions do not distinguish between 2D and 3D, con- 
cerning, e.g., the tendency toward velocity and magnetic 
field alignment or the spectral properties. Thus it is not 
surprising, that numerical studies of MHD turbulence 
have been mostly concentrated on two-dimensional sim- 
ulations, where high Reynolds numbers can be reached 
much more readily, see e.g., [pj- While 2D simula- 
tions are now being performed with up to N 2 = 4096 2 
modes (or, more accurately, collocation points) stud- 
ies of 3D MHD homogeneous turbulence have to date 
been restricted to relatively low Reynolds numbers using 
typically TV 3 = 64 3 modes, e.g., 0, fjj, precluding an 
inertial range scaling behavior. Also in Ref. 0, where a 
somewhat higher Reynolds number could be reached by 



using 180 3 modes, attention was focussed primarily on 
the process of turbulence generation from smooth initial 
conditions and the properties of the prominent spatial 
structures, current and vorticity sheets. 

In this Letter we present results of a numerical study 
of freely decaying 3D MHD turbulence with spatial res- 
olution up to 512 3 modes. We discuss the decay laws of 
the integral quantities, in particular the energy E and 
the ratio of kinetic and magnetic energies T — E v /E M , 
and their dependence on the quasi-constant value of H. 
The energy decay is found to follow a simple law, which 
is determined by T(t) and H. While most previous stud- 
ies have been restricted to the case of negligible magnetic 
helicity H ~ 0, we focus attention on the properties of 
the turbulence for finite H, which is more typical for 
naturally existing MHD turbulence occuring mostly in 
rotating systems. We find that for finite H the energy 
decays significantly more slowly than for H ~ 0. This 
behavior is primarily caused by the rapid decrease of the 
energy ratio T, which has the same decay time as the 
energy. 

The 3D incompressible MHD equations, written in the 
usual units, 



3 t B-Vx(vxB) = ^(-lf^V^B 



(1) 



<9 t w — V x (v x w) — V x (j x B) = //„(- 



w = Vxv, j = V x B, 



(2) 



are solved in a cubic box of size 27r with periodic bound- 
ary conditions. The numerical method is a pseudo- 
spectral scheme with spherical mode truncation as con- 
veniently used in 3D turbulence simulations (instead of 
full dealiasing by the 2/3 rule chosen in most 2D simula- 
tions). Initial conditions are 



B k =ac~ k2 / k °- i 



Vk = be 



— k /k o —i0 k 



(3) 



which are characterized by random phases cvk, Pk and 
satisfy the conditions k • Bk = k • Vk = as well as 
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E = 1 and T = 1. Further restrictions on Bk and Vk 
arise by requiring specific values of H and K, respec- 
tively. The wavenumber ko, the location of the maxi- 
mum of the initial energy spectrum, is chosen as fco = 4, 
which allows the inverse cascade of to develop freely 
during the simulation time of 10-20 eddy turnover times. 
This implies a certain loss of inertial range, i.e., a reduc- 
tion in Reynolds number, but the sacrifice is unavoid- 
able in the presence of inverse cascade dynamics. Choos- 
ing ko ~ 1 would lead to magnetic condensation in the 
lowest-fc state, which would affect the entire turbulence 
dynamics. We have used both normal diffusion v = 1 
and hyperdiffusion v — 2. Apart from the fact that in- 
ertial ranges are wider and H is much better conserved 
for v = 2 than for v = 1, no essential differences are 
found between the two cases. The generalized magnetic 
Prandtl number r\ v j ' fi^ has been set equal to unity. Table 
| lists the most important parameters of the simulation 
runs. 

The energy decay law is a characteristic property of 
a turbulent system. In hydrodynamic turbulence the 
decay rate depends on the energy spectrum at small 
k. Assuming time invariance of the Loitsianskii integral 
£ = f^dl l 4 (vi(x+l)vi{x))the energy has been predicted 
to follow the similarity law E ~ <~ 10 / 7 0. The invari- 
ance of C has, however, been questioned, see e.g., Q. 
Both closure theory [ fl2"| and low-Reynolds number simu- 
lations H yield a significantly slower decrease, E ~ t~ l . 
Experimental measurements of the energy decay law t~ n 
are rather difficult and do not give a uniform picture, n 
ranging between 1.3 Q and 2 p[ . 

The invariance of the Loitsianskii integral has recently 
also been postulated for MHD turbulence ||], where 
£mhd is defined in analogy to C in terms of the lon- 
gitudinal correlation function (z ; ± (.x + Z)z ; ± (x)} of the El- 



saesser fields z 



± _ 



v±B. Since 



E, this assumption 



gives £mhd ~ L°E = const, where L is the macroscopic 
scale length of the turbulence. In addition the expres- 
sion for the energy transfer dE/dt = — e ~ —z 4 /LBq 
was used, which formally accounts for the Alfven effect 
@, 0. These relations give (dE/dt)B / E ll / V ° = const 
and hence E ~ <~ 5 / 6 , treating B as constant. One 
may, however, argue that the Alfven effect is only im- 
portant on small scales I <C L, while on the scale L 
of the energy-containing eddies Bq is not constant but 
B ~ E 1 ! 2 (except for the case that B a is an exter- 
nal field, which would, however, make the turbulence 
strongly anisotropic), hence e ~ E 3 / 2 /L, which would 
give the same result n = 10/7 as predicted for hydrody- 
namic turbulence. Low-resolution numerical simulations 
Q indicate n ~ 1, which is also found in recent simula- 
tions of compressible MHD turbulence |l8) . 

For finite magnetic helicity H provides a constant dur- 
ing energy decay, which for high Reynolds number is 
more robust than the questionable invariance of the Loit- 



sianskii integral. It is true that in contrast to the 2D 
case, where E and are tightly coupled, such that 
E M ^ implies ^ 0, in 3D a state with H = 
and finite magnetic energy is possible. But this is only 
a special and not typical case, since in nature magnetic 
turbulence usually occurs in rotating systems, which give 
rise to finite magnetic helicity. 

If the process of turbulence decay is self-similar, which 
also implies that the energy ratio T remains constant, the 
energy decay law follows from a simple argument |l5[ ]. 
With the scale length L = E 3 ! 2 jc, the dominant scale of 
the energy-containing eddies, we have 

H ~ E M L ~ EL, (4) 



EL, 

since owing to the assumed self-similarity E AI 
E. Inserting L gives 



E v 



dE 
~dt 



£5/2 



(5) 



which has the similarity solution E ~ i~ 2 / 3 . In Fig. 
the ratio E b / 2 /(eH) is plotted for the runs from Table 
with H ^ and small initial correlation po- The figure 
shows that this quantity is not constant, but increases in 
time. Moreover, there is a significant scatter of the differ- 
ent curves. Integration yields a slower asymptotic energy 
decay than predicted n ~ 0.5 — 0.55. (The log-log repre- 
sentation of E(t), often given in the literature to make a 
power law behavior visible, is misleading, since the major 
part of such a curve refers to the transition period of tur- 
bulence generation. The solution (t — t tf )~ n approaches 
the power law t~ n only asymptotically for t^>t*, where 

is not accurately known. We therefore prefer to plot 
the decay law in the primary differential form.) 

We can attribute this discrepancy to the fact that the 
turbulence does not decay in a fully self-similar way. 
Indeed the energy ratio T is found to decrease rapidly, 
in contrast to the 2D case, where T decays much more 
slowly, typically logarithmically g], |J. (The ratio of 
viscous and resistive dissipation e^/e 71 , however, remains 
constant just as in the 2D case |J, which simply reflects 
the basic property, that dissipation takes place in current 
sheets and that these are also vorticity sheets, i.e., the 
location of viscous dissipation.) Let us incorporate the 
dynamic change of T in the theory of the energy decay. 
Assuming that the most important nonlinearities arise 
from the v • V contributions in the MHD equations, Eq. 
(0) is replaced by 



(E v y/ 2 - = - 

L (1 



pl/2 £5/2 

+ r) 3 / 2_ H~ : 



(6) 



using the relation (0). Figure || shows that 
(,E 5 / 2 /e#)L 1/2 /(l + L) 3 / 2 is indeed nearly constant for 
t > 2, when turbulence is fully developed, and the scat- 
ter in Fig. |l| is strongly reduced. Hence relation (|6|) is 
generally valid for finite magnetic helicity. It is also in- 
dependent of the magnitude of the dissipation coefficients 



2 



and character of the dissipation (y = 1 or 2), as long as 
H is well conserved. 

Also the time evolution of the energy ratio T exhibits 
a uniform behavior which is demonstrated in Fig. ||. The 
slight shift of the uppermost curve corresponding to the 
smallest value of H (run 4), is due to the smaller drop 
of r during the very first phase of turbulence generation 
t < 0.5 not included in the figure. Moreover, we find that 
r(i) is proportional to E(t), T ~ cE/H, c = 0.1 - 0.15, 
as seen in Fig. ||, where T/(E/H) is plotted. Inserting 
this result in Eq. (Q) we obtain the differential equation 
for E, which in the asymptotic limit r< 1 becomes 

with the similarity solution E ~ t~ - 5 . For finite T 
the theory predicts a somewhat steeper decay flatten- 
ing asymptotically to t~ ' 5 as F becomes small, which 
is exactly the behavior of E(t) observed in the simula- 
tions. (Note, that if E(t) is plotted on the traditional 
log- log scale, which overrates the transition period t ~ 1, 
a steeper decay would be suggested.) The relation V oc E 
now gives also the similarity law for the kinetic energy 
E v ~t-i. 

This theory does not apply to the special case H = 0. 
Here we find indeed a different decay law, E ~ t~ l from 
run 3, which is consistent with previous simulations at 
lower Reynolds numbers Q and with the prediction in 
Ref. 0. The transition to the slower decay for finite H 
occurs at relatively small values, 0.1-0.2 of the maximum 
possible value. 

We have also studied the effect of an initial velocity and 
magnetic field alignment po = K/E. For small po < 0.1 
the alignment, after increasing initially, tends to saturate 
at some small value, which is due to the fact that K is less 
well conserved than H. For higher po > 0.3 (runs 9 and 
10 in Table I) the alignment becomes very strong, which 
as expected slows down the energy decay drastically. 

In conclusion we have presented a new phenomenol- 
ogy of the energy decay in 3D incompressible MHD tur- 
bulence, which agrees very well with direct numerical 
simulations at relatively high Reynolds numbers. Wc 
consider in particular the case of finite magnetic he- 
licity H, which is typical for naturally occuring mag- 
netic turbulence. The energy decay is governed by the 
conservation of H and the time evolution of the en- 
ergy ratio T = E v / E M . We find that the relation 
(E 5 ^/eH)T 1 / 2 /(l + T) 3 / 2 ~ const is satisfied for most 
H- values and is independent of the magnitude of the dis- 
sipation coefficients and the order of the diffusion oper- 
ator, provided the Reynolds number is sufficiently high 
such that H is well conserved. The kinetic energy is 
found to decrease more rapidly than the magnetic one, 
in contrast to the behavior in 2D, in particular we find 
r oc E. This proportionality leads to a simple energy 
decay law, —dE/dt ~ E 3 , or E ~ t~ - 5 . We also obtain 



the similarity law for the kinetic energy E v ~ t _1 . For 
the special case H = the energy decays more rapidly, 
E ~ t , which agrees with previous simulations at lower 
Reynolds numbers. The transition to the finite- H behav- 
ior occurs at relatively small values of H. 



Results concerning the spatial scaling properties of 3D 
MHD turbulence will be published in a subsequent paper. 
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in time indicates an energy decrease slower than t 2 ^ 3 , typi- 
cally t" ' 5 . 
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FIG. 2. Energy decay law 

in differential form (E 5/2 /eH )F 1/2 /(l + T) 3/2 ' for the same 
runs as in Fig. |l|. The lowest curve, which falls somewhat 
outside the main curve bundle, corresponds to the run with 
the smallest Reynolds number (run 1), where conservation of 
H is least good. 
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FIG. 4. T/{E/H) for the same runs as in Fig. [j] demon- 
strating the proportionality T ex. E. 



TABLE I. Summary of the simulation runs. The value of 
H — 0.28 corresponds to the maximum value for the given 
spectrum (^), H < E/ko- 



run No 


N 


V 


r\u 


H 


pa 


tmax 


1 


256 


1 


10~ 3 


0.19 


0.04 


18.5 


2 


512 


1 


3 x 10~ 4 


0.19 


0.04 


10 


3 


256 


2 


10~ 6 





0.05 


20 


4 


256 


2 


10~ 6 


0.11 


0.05 


10 


5 


256 


2 


1Q -6 


0.19 


0.04 


20 


6 


512 


2 


3 x 10~ 8 


0.19 


0.04 


10 


7 


256 


2 


10~ 6 


0.25 


0.04 


10 


8 


256 


2 


1Q -6 


0.28 


0.03 


10 


9 


256 


2 


10" 6 


0.19 


0.38 


10 


10 


256 


2 


10~ 6 


0.19 


0.71 


10 
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